Theory of robust subrecoil cooling by STIRAP 
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We demonstrate one-dimensional robust Raman cooling in a three-level A-type atom, where the 
velocity-selective transfer is carried out by a STIRAP pulse. In contrast to the standard Raman 
method, the excitation profile is insensitive to variations in the pulse duration, whereas its position 
and width are nevertheless under control. Two different cooling variants are examined and we show 
that subrecoil temperatures are attainable in both cases. 
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I. INTRODUCTION 

The physics of cold atoms is an important ingredient in 
the emerging field of quantum technology, either provid- 
ing elements for demanding applications or a transparent 
testing ground for ideas and methods |TJ-[3] . This is made 
possible by sophisticated control of both translational 
and internal degrees of freedom for neutral atoms with 
electromagnetic fields, demonstrated especially by laser 
cooling and trapping methods [3HS]. One of the most 
notable successes has been the reaching of quantum de- 
generacy, i.e., Bose-Einstein condensation with bosonic 
atoms, and Cooper pairing with fermionic atoms [SHI]. 
Here the key tools have been evaporative cooling and 
magnetic trapping, as reaching the ultralow tempera- 
tures and high densities with laser light is usually ham- 
pered by heating and loss of atoms due to light-assisted 
atomic collisions [H][I2] or reabsorption of scattered pho- 
tons [T31.IT4]. 

Evaporative cooling, however, is rather wasteful of 
atoms and requires fast thermalizing collisions, which 
limits its practicality in some cases, such as single-species 
fermionic gases [15j [16], atomic clocks [10], or small 
samples of cold atoms [5] , and especially transversal cool- 
ing of atomic beams [TTJ [TS] . Very recently, for instance, 
fermionic atoms have been efficiently laser cooled using 
narrow optical transitions, which lead to low Doppler 
temperatures [131 [IB]. The same approach has been 
used earlier with alkaline-earth atoms, for the purpose 
of building atomic clocks for optical wavelengths [8, 10J. 
In general, state-insensitive traps are needed for quantum 
state engineering and precision metrology [8Tll0|. Simi- 
larly, one needs to develop alternatives also for evapora- 
tive cooling. A simple possibility is to develop further the 
concept of Raman cooling; one of the aspects that can be 
improved is the robustness in respect to the laser pulse 
parameters. One of the success stories in robust quantum 
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control of internal state dynamics is the use of adiabatic 
following of eigenstates, i.e., stimulated rapid adiabatic 
passage (STIRAP) [19 . As known e.g. with molecular 
dynamics [20 , 21] , it can also be used to control motional 
degrees of freedom. 

Raman cooling [22] is one the most efficient non- 
evaporative techniques for optical cooling of atoms be- 
low the one-photon recoil limit. The cooling efficiency 
is associated with population transfer which is typically 
generated by a 7r-pulse that has a Blackman [22H24] or 
a square [25] HBJ envelope. Although such a 7r-pulse op- 
erates efficiently on atoms corresponding to the resonant 
velocity group, the process requires extreme control of 
the pulse amplitude and duration. Even though some 
optimization can be obtained for the method |27j , an ap- 
proach based on adiabatic passage is of definite interest 
due to its robustness regarding pulse amplitudes and du- 
rations. However, the use of adiabatic passage in Raman 
cooling is not possible without having an exact under- 
standing of its velocity-selective properties. In this paper 
we present a careful analysis of the Raman cooling via 
STIRAP pulses, and show that subrecoil temperatures 
are attainable. 

The possibility to use adiabatic passage for optical 
cooling was demonstrated by Korsunsky [21], being ap- 
plied to VSCPT [22] cooling. Raman cooling is a much 
greater challenge which requires a careful analysis of the 
velocity selection process. Another difficulty consists in 
a reasonable choice of adiabatic passage which would 
demonstrate a benefit from its usage in Raman cool- 
ing. For instance, Raman cooling employing rapid adi- 
abatic passage (RAP) was presented in papers [5UH52] . 
where the amplitude and the frequency of Raman pulses 
are changed in a controlled way, producing a velocity- 
selective transfer of the atomic population. The extent 
of the frequency-chirp determines the range of excitation 
by RAP and makes it possible to transfer even the wing 
of the velocity distribution. However, RAP does not give 
an appreciable advantage over ordinary Raman cooling, 
because the excitation profile depends on the amount 
of chirping, the pulse area and its envelope. The the- 
ory presented in this paper shows that the utilization of 
stimulated Raman adiabatic passage (STIRAP) [33J [33] 



2 



in contrast leads to a crucial growth of cooling robust- 
ness. In fact, once the adiabaticity criterion is fulfilled, 
the velocity selection in a case of large upper-level detun- 
ing depends only on the pulse amplitudes, allowing wide 
variations in the pulse duration. 

This paper is organized as follows. After describing 
the basic atom-pulse setup in Sec.|TTJ we discuss the adia- 
baticity criterion in Sec.|III| Assuming that this criterion 



is satisfied, we present in Sec. IV the conditions of the 
adiabatic conversion for internal atomic states required 
for cooling. In Sec.JV] we discuss the velocity selection 
process and in Sec. |VI| the elementary cooling cycle. As 
the full cooling process requires a series of these cycles, 
and the adiabaticity criterion should be fulfilled for each 
cycle, we formulate in Sec. |VII| two possible approaches: 
a) using STIRAP pulses of equal amplitude, but with du- 
rations and amplitudes that vary during the pulse cycle 
series (in analogy with the standard Raman method, see 
e.g. Ref. [57]), and b) using STIRAP pulses of substan- 
tially different amplitudes, but not changing their dura- 
tion and amplitude during the pulse cycle series. In the 
latter case, the amplitudes are adjusted to transfer the 
wing of the velocity distribution together with the possi- 
bility to cool the atomic ensemble below the recoil limit. 
The cooling results for a series of cycles for these two 
different transfer types are presented and discussed. The 
paper is concluded by the summary given in Sec. |VIII| 



II. THE ATOMIC HAMILTONIAN 

Figure [lja) illustrates a three-level A-type atom that 
travels along direction Oz with velocity v, being origi- 
nally prepared in ground state |1). A pump laser prop- 
agating along the Oz axis couples the atomic transition 
|1) <-» 1 2), whereas the contra-propagating Stokes laser 
couples the transition |2) <H- |3). The electric field of 
laser pulses is written as 



E(r,t) = -E P e 



ikz—iojpt 



:E s e 



-ikz — iujst 



where u)p, u>s ar e the corresponding laser frequencies, 
and k is the wave number for both light beams (the rele- 
vant frequency quantities are in fact the atom-field detun- 
ings, which are small compared to the actual frequencies, 
and thus we can assume equal values of k). The pump 
and Stokes pulses produce a velocity-selective STIRAP 
from the |1) to |3) ground state. 

The total Hamiltonian for the atom-light system con- 
sists of the kinetic-energy term P 2 /2M, the Hamiltonian 
of the A-type atom and the interaction Hamiltonian V: 



H = 



2M 



The laser-atom coupling V is the sum of couplings with 
both the pump and Stokes pulse, and in the rotating- 




Rabi 

frequencies 




At 

(b) 



FIG. 1. (a) A three-level A-type atom interacting with a pair 
of contra-propagating laser beams. The pump laser links the 
initial state |1) with an intermediate state |2), whereas the 
Stokes laser links the intermediate state with the final state 
j3). (b) Both the pump and Stokes pulse have a Gaussian 
profile in time, overlapping during period At. 



wave approximation (RWA) is given by 

^21-Ep| oWlU jfe z _i Wp< 



V P = -^^\2)(l\e« 



|2)<3|< 



-ikz — iujst 



H.c. 



where d%i, (I23 are the corresponding electric dipole mo- 
ments; H.c. is the Hermitian conjugate. The Rabi fre- 
quencies 

n P = — — , n s = — — 

are considered to be real-valued and positive, and thus 
the coupling V is written as 



V 



/l ^|2)(l|e lfe2 -^ t + ^|2)(3|e-^-^ st + H.c. 



(1) 



Considering z in Eq. (JlJ as an operator acting on the 
external degrees of freedom of the atom, we apply in 
Eq. ([TJ the relationship 

e ±ik* = J2\ P ){pThk\. 

p 

The obtained laser-atom coupling 



v = ^(^\2, P )(i,p-nk\ ( 



2,p)(3,p + hk\e- lLUst + H.c. (2) 



shows that three states coupled by the pump and Stokes 
pulses form a momentum family 

F(p) = {\i,p-hk),\2, P ),\3,p + hk)}. 
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As long as spontaneous emission is not taken into ac- 
count, a single family J-(p) of a specific momentum p 
can be only taken into account, so the sum in Eq. ([2| is 
avoided. 

In the basis of the three states 



The corresponding (time-dependent) dressed-state eigen- 
values are 



ax 



exp 



\a 2 ) = exp 
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(p - hkf 


h 


2Mh 




(p-hkf 
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2Mh 


E 3 


( P - hkf 
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2Mh 
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\l,p-hk), 
2,p), 



Ap-A, 



the dynamics of an atom starting from state |<ii) with 
the corresponding velocity v = (p — hk)/M is described 
by the atomic Hamiltonian 



H 







1 



n P (t) 



o 



i 



\n P {t) Jt-u R -A P - 

1^ / n 2fcn 
-fist) — +A S 

2 6W M 



A, 



(3) 



where Ap = tup — W21, A5 = — W23 are the pulse 
detunings; top — hk 2 /2M is the recoil frequency. 

The pump and Stokes pulses evolve in time as Gaus- 
sian profiles shown in Fig. [ljb), being arranged in the 
counterintuitive sequence with t$ <tp: 

Q P (t) = JW- (t -^ )2/2Tj % n s (t) = n so e-^ 2 /* T s, 

where 2Tp, 2T$ are pulse widths. These lasers stimulate 
a velocity-selective STIRAP from the \a\) state to the 
|d3) state. The corresponding resonant velocity group 
is considered in Sec. whereas Sec. |IV| and Sec. [V| are 
devoted to an arbitrary velocity group in the case of large 
enough detunings Ap and A5. 



III. TWO-PHOTON RESONANCE 

Two-photon resonance is characterized by momentum 
p = M(Ap — As) /2k, for which the detuning from state 
1 3) in representation ^ equals zero. Thus an atom pre- 
pared in state |1) falls into the resonance, if it starts with 
resonant velocity 



A P - A s 
v = 2jfe Wrec ' 



(4) 



where v xec is the recoil velocity. Under this condition, 
the basis of time-dependent eigenstates of the Hamilto- 
nian ([3| is given by (see Ref. [34]) 

\a + ) = sin 9 sin $|ai) + cos^ja^) + cos 9 sin $| 03), 



\a ) = cos9|ai) — sin9|a 3 ), 

\a~) = sin9cos$|ai) — sin$|a2) + cos 9 cos $| 03). 



(5) 



a 2 p + n 2 P (t) + n 2 s (t) - A P , uj° = o, 



A 2 P + n 2 p (t) + n 2 (t) + Aj 



where Ap = Ap + cop — kp/M . The angle $ is a function 
of the Rabi frequencies and detunings |35j : 



|3,p + hk), 



tan $ = 



y/n 2 p (t) + n 2 s (t) 



Ap 



n 2 p (t) + n 2 (t) - a p 



(6) 



whereas the mixing angle 9 depends only on Rabi fre- 
quencies: 



tan 9 



(7) 



In the case of large detuning Ap (fip(t), Qs'(i) <C 



|Ap|), the basis of eigenstates {\a^ 



\a )} is re- 



duced to the basis of states {|C), |NC), \a 2 )}- Both the 
coupled |C) and non-coupled |NC) state is a combination 
of ground states: 



|C) = sin9|ai) + cos9|a3), 
|NC) =cos9|ai) -sin9|a 3 ), 



(8) 



with the corresponding eigenvalues 



n 2 P (t) + n 2 {t) 

4Ap" : 



wnc = 0. 



Here, we take into account that |Ap| 3> cop, kp/M . 

The atomic population is contained in states |C) and 
|NC), the Hamiltonian matrix element for nonadiabatic 
coupling between these states is given by (C|^|NC) [55] , 
The "local" adiabaticity constraint reads that this matrix 
element should be small compared to the field-induced 
energy splitting \u c ~ wnc|, i-e., 



|(C|-|NC)| « I^c-^ncI 



191 « 



n 2 P (t) + n 2 (t) 



Taking a time average of the left-hand side, (9 a ^) = 
7t/2At, where At is the period during which the pulses 
overlap, one gets a convenient "global" adiabaticity cri- 
terion 



n 2 p {t) + n 2 {t) 



At > 1. 



(9) 



IV. 



THE EIGENSTATES OF THE EFFECTIVE 
HAMILTONIAN 



In the case of large detuning Ap (flp(t),ilg(t) <C 
|Ap|), adiabatic transfer can be considered at arbitrary 
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velocity. Then the upper state \a 2 ) is almost unpopulated 
during STIRAP and can be adiabatically eliminated. 
Hence, substituting the probability function \^f) into the 
Schrodinger equation, one may assume that (02! 4i\^) ~ 
0. Contributions into state \a 2 ) are derived from the 
Hamiltonian ^ and, in the case of |Ap| ^ wjj, kp/M, 
are given by 



(a 2 m 



n P (t) 

2A P 



ail*) 



2A P 



(as|*). (10) 



Approximation ( 10 ) requires that the condition 



|( a2 |||*)|«|A P ||(a 2 |*)|, 



(11) 



is fulfilled. By substituting Eq. ( 10 ) into the left-hand 
side of inequality ( 11 ), one gets the necessary constraints 



|A s -Ap|,T- 1 «|A P |. 



Using Eq. (10), STIRAP can be described in the basis 
of states {|ai)7]a3)}, and the effective Hamiltonian of this 
two-level system is written as 



-2S (t) n cS (t) 

Q eS (t) 2(5 eS {t) - S (t)) 



(12) 



The effective detunings and the Rabi frequency are 

n P (i)n s (t) 



s (t) 



nf,(t) 



4A, 



n cS (t) 



2A, 



5 e g(t) = AS- 



Q|(t) - n%(t) 

4A~p" : 



(13) 



where AS = A s - A P + 2kp/M = 2k(v - v ); v is the 
original velocity of an atom along axis Oz. The effective 
Hamiltonian (12) has the following eigenstates: 



\a + ) = sin <d\ai) + cos 9|a 3 ) 
|a_) = cos 9|ai) — sin 9|a 3 ) 



tan 9 




with the corresponding cigenfrcqucncies 
w± = -S 



s e s fw 
2 ~2~ 



1 



cir 



The approach of eigenstates \a + ) and |a_) covers the 
case of two-photon resonance Q as a specific case with 
A5 = 0. In this case, the mixing angle 9 in Eq. ( 14) coin- 
cides with that in Eq. frfh, and eigenstates |a + ), a_) turn 
into states |C) and |NC), respectively. During STIRAP, 
the angle 9 varies from to tt/2, and the non-coupled 
state |NC) evolves from the |ai) to \ag) state, involving 
an adiabatic transfer of population. Consequently, the 
eigenvalue wnc changes from the frequency of state \ai) 
to that of state (03). 

General case of arbitrary A<5 can be considered in the 
same way where the role of state |NC) is played by the 
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FIG. 2. The eigenfrequencies of pure atomic states |ai) and 
I03} (dashed lines) cross at point 5 e g = 0. The avoided level 
crossing for the eigenstates (solid lines) leads to a correspond- 
ing change between the eigenstates and pure states. 



eigenstate |o_). STIRAP occurs together with a change 
in the eigenfrequencies ui± , whose dependence on detun- 
ing 8 c g is shown in Fig. [2} While the laser pulses overlap, 
<5 e ff can vary from negative to positive value or vice versa, 
in dependence on the sign of Ap. The switching point 
between states |ai) and \as) corresponds to a crossing be- 
tween the frequencies of these pure states, which appears 
at (S e ff = 0. 

The time evolution of the eigenfrequencies ui± is shown 
in Fig. [3] where three possible cases against the detuning 
<5 c ff are presented: (a) 5 c g > or (b) S e g < all the 
time while laser pulses overlap; (c) S e g = at a specific 
time during this period. Both eigenstates \a+) and |a_) 
in Figs. |3][a) and [3jb) return to the initial pure state 
in the end of STIRAP, suppressing adiabatic transfer of 
atoms. A crossing between the frequencies of states |oi) 
and \a,3) appears in Fig. [3jc), where atoms that adia- 
batically follow the eigenstate |a_) are transferred from 
the |ai) to ^3) state. The corresponding velocity group 
derived from the explicit form ( 13 ) of 5 c g is given by 



"50 
~8fc|A P 

~8fc|A P 



< V — Vo < 



< V — V < 



n 2 

"P0 

8k\A P 

O 2 
"so 

8k\A P 



if Ap > 0, 



if Ap < 0. 



(15) 



If Rabi frequency f2 e g(t), detunings 6 c g(t) and So(t) 
vary sufficiently slowly so that the nonadiabatic mixing 
between states |a+) and |a_) is negligible, population re- 
mains in state |o_) during the whole STIRAP pulse. If 
this condition holds for all velocity groups, the velocity 
range (15) corresponds to population transfer from state 
|ai) to 103). The transfer is accompanied by the probabil- 
ity near unity for all atoms within the range (15), is not 
affecting outside atoms. However, because the transfer- 
should be adiabatic for huge velocities, required pulse du- 
ration is so long that an acceptably efficient cooling can 
not be achieved. In contrast, we suggest the adiabaticity 
criterion to be satisfied only for atoms at the two-photon 
resonance Q, which in turn allows shorter durations of 
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is valid regardless of whether the transfer is adiabatic or 
not. If atoms follow the eigenstate |a_) adiabatically, 
both population in state |a 3 ) and the mixing angle 
has the similar behaviour, because the population varies 
as sin 2 0. We assume that has the largest rate for 
resonant atoms as well, at least for quite large At. This 
fact relative to an average value (0 O u) gives an inequality 

(@av) < tt/2At. 




las)-**. 



Ar 



(b) 



l«a) 



(c) 



In cases shown in Figs.[3]Ja) and[3]Jb), in spite of the fact 
that each atomic population at the beginning and the 
end of laser pulse takes the same value, can change 
significantly in range from to 7r/4. Hence (Q a v) can be 
of the same order of magnitude as At -1 . 

After the substitution of S e s(t), fi e ff(i), to the inequal- 



ity (16 1, the adiabaticity criterion reads 



At 



'fip(t) 



4A, 



AS AS 



2A, 



(17) 



The two-photon resonance corresponds to Ai5 = 0, 



when the inequality (17 1 coincides with the adiabaticity 
criterion ([9]). As seen from Eq. (17), once the adiabatic- 
ity criterion ^ is fulfilled, atoms following state |a_) 
adiabatically are given by 



FIG. 3. The time evolution of the dressed-state frequencies 
(solid line) on the pulse-overlapping region At. The dashed 
lines represent the frequencies of the pure states. In depen- 
dence on AS, one gets the following cases: (a) <5 c ff > 0, (b) 
5 e ff < during this period; (c) S c ff — at a specific time, (c) 
represents adiabatic transfer from the |ai) state to |as) state. 



the STIRAP pulse. The corresponding resonant velocity 
range for STIRAP is considered next in the next section. 



AS AS 



n|(t) 



n 2 P (t) 



2A f 



> 0. 



Hence, in addition to the resonance-velocity group, adi- 
abatic transfer includes the following velocity groups: 



v < 



v < 



n 2 

"so 
'4fe|A P | 

O 2 

"P0 

"4fc|A P | 



or v 



v > 



n 2 

"po 
4fc|A P 

o 2 
"so 

4fc|A P 



if A P > 0, 



if A P < 0. 



V. VELOCITY SELECTION OF STIRAP PULSE 

In the case of large detuning Ap, the state vector 
follows the dressed state |a_) adiabatically throughout 
the interaction if a condition 



\{a+\j t \a-)\ « |«- 



is fulfilled during the STIRAP pulse. After the substitu- 
tion of states \a+) 1 |a-) and the corresponding eigenfre- 
quencies from Eq. (14 1, this constraint is written as 



iei« 



"e 2 ff 



(16) 



Next we consider the "global" adiabaticity criterion, re- 



placing Q in Eq. ( 16 1 by an appropriate time average 

The rate of population transfer from state |ai) to I03) is 
the largest for atoms in the two-photon resonance, which 



However, in this case, equals at the end of STIRAP 
pulse, and atoms return to state |ai). That is why only 
those atoms can be transferred to state \as), for which 



O 2 
"so 

~4fc|Ap 

"po 
4fc|Ap 



< v — v < 

< V — v < 



O 2 
"po 

Ak\A P 

n 2 

"so 
4k\A P 



if Ap > 0, 



if Ap < 0. 



(18) 



The velocity range ( 18 ) is only defined by the two- 



photon Rabi frequencies, being twice broader than that 
by Eq. (15 1. The transfer probability reaches unity at 



the two-photon resonance and then declines to zero as 
the boundary of range ( 18 ) approaches. Possible veloc- 



ity envelopes will be discussed in Sec. |VII| as applied for 
Raman cooling. 

The STIRAP pulse transports population from the 1 1} 
to 1 3) state without an occupation of the immediate state 
1 2) for the resonant-group atoms only. For the rest of 
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atoms, the transfer is accompanied by populating the |2) 
state, which in turn gives rise to a spontaneous decay 
from this state to ground states of the A-type atom. 

A contribution from spontaneous decay can be esti- 
mated with help of the density operator <r, whose matrix 
elements are given by 



<Tij(p) = (ai\<r\a,j), i,j = 1,2,3. 



(19) 



Portions of atoms in state \a 2 ) and those of them leaving 
state \a 2 ) due to spontaneous decay, 022 (p) and <x sp (p), 
are related by relationship 



d 
dt 



<j sp (p)=T J CT 2 2(p-ftfcn-e z )$(n)d 2 n, 



where T is the rate of spontaneous emission, d 2 n is the in- 
finitesimal solid angle in direction n of spontaneous emis- 
sion. Because state \a 2 ) is only populated during period 
At, while the pump and Stokes pulses overlap, <r sp (p) is 
given by 



<?s P (p) 



ry J o- 22 (p- hkri ■ e z )<S>(n)d 2 ndt. (20) 



The inequality ( 10 ) imposes a constraint on a fraction of 



atoms m state |a 2 ) 



0"22(PJ < 



2A P 



o-ii (p) 



2kJ aM ^ 2A| • 



Substituting 022 (p) into Eq. ( 20 ) and taking into account 
that J $(n)(i 2 n = 1, one gets a contribution of the spon- 
taneous decay 

MP) < r ^2 At. 

Atoms spontaneously decaying from state |2) are ne- 
glected from the consideration if a sp (p) <C 1. This con- 
dition in combination with the adiabaticity criterion Q 
gives the constraint 

> k 1 Ar > 1, 



r |A P | 

which can be reached by increasing detuning Ap. 

VI. ELEMENTARY COOLING CYCLE 

The first step of a cooling cycle consists in transferring 
atoms from state |1) to |3) due to STIRAP pulse. Once 
the adiabaticity criterion |9]) is fulfilled, atoms at the 
resonant velocity vq Q are transferred with efficiency 
near unity. For laser configuration depicted in Fig. [T] 
atoms after the transfer get a momentum gain of 2hk 
along Oz, so one should follow with a condition i> < 
0. Because the spread of the resonant-velocity group is 
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FIG. 4. (a) The excitation probability (colored online) of 
atoms in state |1) against velocity v by STIRAP pulses of 
different coupling parameters according to set (22 1, (b) The 
velocity spread Av decreases from the initial 3u rC c to 0.1« re c 
of cooled ensemble. 



defined by Eq. ( 18 ), atoms at zero velocity are suppressed 



from the transfer when 

if A P > 0, 



< - 



"P0 



4fe|Ai 



"0 



< 



O 2 
"so 

~4fc|A P 



if A P < 0. 
(21) 



The inequality ( 21 ) only imposes a constraint on maximal 



magnitude of Rabi frequencies and thereby of the veloc- 
ity spread of the resonant atoms. By decreasing Rabi 
frequency Qso or f2po, the magnitude of the resonant ve- 
locity can be chosen substantially smaller than the recoil 
velocity. As a result, STIRAP can be efficiently used as 
the first step of an elementary cooling cycle. But such 
a transfer does not need the exact holding of the pulse 
duration, as that does in the standard Raman cooling 

In the second step, optical pumping excites atoms from 
1 3) to 1 2) state, changing the momentum of an atom from 
p' to p' — hk. Then the atom spontaneously decays in 
the initial state |1), and its momentum becomes p = 
p' ' — Hk— Ap where Ap is the projection of a spontaneously 
emitted photon; |Ap| < hk. The resulting population in 
state |1) is written as 

(l,p|cr'|l,p) = (l,p\a\l,p) + (3,p+hk+Ap\cr\3,p+hk+Ap), 



which with help of Eq. ( 19 1 takes the form 



°ii(p) = ct ii(p) + °33(P -Hk + Ap). 
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FIG. 5. (a) The excitation probability of atoms in state |1) 
against velocity v by STIRAP pulse ( 23a I with the maximum 
at the resonant velocity Vq = — 0.25iw. (b) The velocity 
spread Av decreases from the initial 3w roc to 0.1« rcc of cooled 
ensemble. 



After the cooling cycle the number of atoms at zero ve- 
locity increases together with the cooling of the atomic 
ensemble. 



VII. COMPARISON OF TWO COOLING TYPES 

In this section, we compare two different variants of 
cooling by STIRAP. The first one uses the similar laser 
beams with Clpo — Qso, so the velocity profile of the 
pulse has a symmetric shape. A condition when atoms at 



zero velocity are not transferred is derived from Eq. (21 1, 
and is written as 



n% = n| = 4fc|u ||Ap|. 

Both the pump and the Stokes pulse has the same pulse 
shape, so 

T P = T s = 0.5(tp~t s ). 

As the resonant velocity v gets closer to zero, the width 
of the excitation profile decreases. To entirely excite the 
wing of the initial profile of spread 3i> roc , we use a set of 



I v \ =4*W 2 * 



0, 



(22) 



for which the probability of transfer from state |1) is 
shown in Fig. [4^a) . Note that as |«o| decreases, the mag- 
nitudes f2p /|Ap| and 57| /|Ap| decrease as well, which 



in turn leads to an increase of Ar in order to fulfil the 
adiabaticity criterion After the set has been applied, 
laser beams are alternated in order to excite the right- 
side profile. Such a cooling method is similar to ordinary 
Raman cooling. 

In the second variant, the complete wing of initial dis- 
tribution except a narrow peak at zero velocity is trans- 
ferred by one STIRAP pulse as shown in Fig. [5ja) . The 
peak width is defined by the resonant velocity Vo which is 
equal to — 0.25w rcc . Notice that such an excitation is not 
attainable for ordinary Raman cooling, and represents a 
new type of a velocity-selective transfer with an essen- 
tially asymmetric profile. The advantage is in utilizing 
a single profile instead of the set ( 22 ) . Rabi frequencies 
are given by 



fip = 4fc|uo||Ap|,fiso=8f2po ifA P >0, (23a) 
n| =4fc|uo||Ap|,Qp =8O S0 ifA P <0. (23b) 



Because a pulse of the largest amplitude should decrease 
faster, pulse widths are given by 

T P = 0.5(t P - t s ), T s = 0.35(t P - t s ) if A P > 0, 
T P = 0.35(ip -t s ),T s = 0.5(t P -t s ) if A P < 0. 

After each STIRAP pulse, laser beams are alternated 
giving the same excitation of the right side of the velocity 
distribution. 

In both cooling variants, the duration T pu i sc of 
STIRAP pulse is defined by the start and end time 



ts - (tj 



''end 



tp + (tp-t S ), 



being equal to 3(<p — ts). Pulse durations in the first 



variant in accordance with the set ( 22 ) are given by 



Tpuisc = 6 • 2 fe rp fc = 0,...,4, 

where tr — Wp 1 is the recoil time. In the second variant, 
each pulse has the same duration equal to 96rp. 

Figures (4^b) and [5jb) show results of both cool- 
ing methods starting from the initial distribution with 
the spread of 3 p rcc . The velocity spread Av = 
(FWHM)/\/8Tri2 has been reduced to nearly a) 0.06u rcc 
in Fig.Qb), and b) 0.04w roc in Fig.[5^b). The correspond- 
ing temperatures of the atomic ensemble have gone down 
to a) 0.004T roc , and a) 0.002T rec , where T rcc is the recoil- 
limit temperature. The duration of optical pumping is 
considered as a negligible value against that of STIRAP 
pulse, which in turn gives total durations of cooling meth- 
ods as 7440rp and 7680rp, respectively. In spite of the 
fact that both durations are approximately equal, the 
number of cooling cycles are rather different: the first 
method contains N = 200 elementary cycles, whereas 
the second one contains N — 80 cycles. For instance, 
being applied for cesium atoms (wp ~ 2ir x 2kHz), both 
cooling variants take a time of 0.6s. 
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VIII. CONCLUSION 



We have described Raman cooling by velocity-selective 
STIRAP in the case of large upper-level detuning. We as- 
sume that the adiabaticity criterion is fulfilled for atoms 
in the two-photon resonance, so these atoms are trans- 
ferred with efficiency that approaches unity. The po- 
sition and width of the excitation profile is insensitive 
to the pulse duration, being in a linear dependence on 
the two-photon Rabi frequencies. This profile can take 
an asymmetric form, exciting the wing of the velocity 
distribution with the exception of a narrow peak near 
zero-velocity group. We have compared such a type of 
STIRAP with that of symmetric profile where a set of 
different coupling parameter utilized for excitation of the 
wing of the velocity distribution. Both methods has well- 
controlled the excitation probability, making the attain- 
able temperatures essentially go below the one-photon 
recoil limit. 

Our simplified treatment does not necessarily take into 
account many aspects of an actual experiment, but it 
shows that a priori the idea of Raman cooling combined 
with STIRAP offers a tool for reaching ultracold tem- 
peratures without the use of evaporative cooling. This 
is an example of robust quantum control of transla- 
tional atomic degrees of freedom. If subjected to trapped 
atoms, ID cooling eventually leads to cooling at higher 
dimensions, although at tight confinement the quantiza- 
tion of motion opens the door for sideband cooling (which 



also relies on a Raman process) [8,37. It can also be ap- 
plied to low-dimensional systems prepared with optical 
lattices [35], such as the elongated quasi-lD cigar-shaped 
clouds in which the lattice provides tight confinement 
in perpendicular direction but atoms are almost free to 
move in axial direction, in which one can then apply Ra- 
man cooling. If using such a system as a wave guide, 
the trapping potential in axial direction is quite weak, 
which limits the use of side-band cooling in axial direc- 
tion. Also, for atomic clocks one needs high numbers of 
atoms and yet low densities to avoid interaction-induced 
frequency shifts, and one solution is again an optical lat- 
tice that slices the sample into non-interacting 2D pan- 
cakes [10]. Although then narrow-band cooling is the 
basic tool for e.g. alkaline earth atoms, one can also con- 
sider cooling at higher-lying atomic energy levels, which 
gives a rich level structure and opens the possibility for 
Raman cooling with tripod configuration. Finally, we 
note that the very recent interest in sub-Dopplcr cooling 
of fermionic isotopes [Ml HO] provides yet another appli- 
cation for non-evaporative cooling methods. 
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